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Abstract 

We introduce the differential and the total decay widths of a resonant (Gamow) 
state decaying into a continuum of stable states. When the resonance has several 
decay modes, we introduce the corresponding partial decay widths and branching 
fractions. In the approximation that the resonance is sharp, the expressions for 
the differential, partial and total decay widths of a resonant state bear a close re¬ 
semblance with the Golden Rule. In such approximation, the branching fractions 
of a resonant state are the same as the standard branching fractions obtained 
by way of the Golden Rule. We also introduce dimensionless decay constants 
along with their associated differential decay constants, and we express experi¬ 
mentally measurable quantities such as the branching fractions and the energy 
distributions of decay events in terms of those dimensionless decay constants. 
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1 Introduction 


The standard Golden Rule allows us to calculate the transition rate (probability of 
transition per unit time) from an energy eigenstate of a quantum system into a con¬ 
tinuum of energy eigenstates. If \E[) is the eigenstate of an unperturbed Hamiltonian 
H 0 , and if such state is coupled to a state \Ef) by a perturbation V, the transition 
probability per unit of time from the initial state | E\) to the final state \E{) is given, 
to first order in the perturbation, by 

f = y \ {E{\V\Ei)\ 2 5(Ei - Ei ), (1.1) 

where (E{\V\E[) is the matrix element of the perturbation V between the final and 
initial states. If the initial state is coupled to a continuum of final states |Ef), and if 
the density of final states (number of states per unit of energy) is p(E{), the transition 
probability per unit of time from the state \E{) to the continuum of final states \Ef) is 
given by 

ifc-f = I d E { p(E { )R^ { = Y \(E { \V\E^\ 2 p(E -^, (1.2) 

where it has been assumed that the matrix element (Ef\V\Ei) is the same for all the 
states in the continuum. The decay width from the initial state | Ei) to the final state 
\E{) is given by 

Ti-rf = hRi^t = 2n\(E { \V\E i )\ 2 6(E { - E f ). (1.3) 

The decay width from the initial state | E\) to a continuum of final states \Ef) is given 
by 

I • = hR^ = 2tt {{EflVlEM 2 p(Ei ). (1.4) 

The standard derivation of Eqs. (Il.lj) - (ll.4j) is the result of first-order perturbation 
theory, and it is valid when the initial and final states are stable. The purpose of the 
present paper is to introduce the analog of Eqs. (Il.lj) - (ll.4j) under the assumption that 
the initial state is described by a resonant (Garnow) state, and that such state decays 
into a continuum of stable, scattering states. 

In Sec. El we introduce the differential and the total decay widths of a resonant 
state that has only one decay mode. In Sec. [3l we derive an expression for such decay 
widths in terms of a truncated Breit-Wigner (Lorentzian) lineshape and the matrix 
element of the interaction. We also show that, when the resonance is sharp and far 
from the energy threshold, the expressions for the resonant decay widths bear a strong 
resemblance with Fermi’s Golden Rule. In Sec. El we apply the results of Sec. El to the 
delta-shell potential. In Sec. El we introduce the partial decay widths of a resonant 
state that has more than one decay mode, we define the branching fractions for each 
decay mode, and we point out that, at least in principle, such branching fractions 
afford a way to falsify the formalism of the present paper. In Sec. El we introduce 
dimensionless partial and total decay constants along with their associated differential 
decay constants, we express the branching fractions in terms of them, and we argue that 
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the differential decay constant corresponds to a measurement of the energy distributions 
of decay events (the invariant mass distributions of particle physics). Section [7] contains 
our conclusions. 


2 Preliminaries 


Let H = Ho + V be a Hamiltonian that produces resonances, where Ho is the free 
Hamiltonian and V is the interaction potential. The continuum spectra of both H and 
Ho will be assumed to be [0, oo), as it often occurs in non-relativistic quantum mechan¬ 
ics. For simplicity, we will assume that the continuum spectra are non-degenerate, and 
that the eigenstates of H 0 can be determined by a single quantum number, the energy 
E. The eigenstates of the free Hamiltonian will be denoted by | E), 

H 0 \E) = E\E). (2.1) 

We will denote the resonant state by |z R ), 

H\z R ) = £ r |zr) , (2-2) 


where z R = E R — iT R /2 is the complex resonant energy of a decaying state. We 
will denote by e~ lHr ^ h \z R ) the time evolution of the resonant state, r being the time 
parameter that appears in the Schrodinger equation. 

When | E) is delta normalized and 1%) is normalized according to Zeldovich’s nor¬ 
malization, the “scalar product” (E\e~ lHr ^ H \z R ) has dimensions of (^/energy) . Thus, 
|(P|e _li?T//ft | 2 : R )| 2 has dimensions of (energy) -1 , and it can be interpreted as a prob¬ 
ability density, i.e., a probability per unit of energy. This is why we will identify 
|(.F|e -1 ' HT ' // *|z R )| 2 with the probability density ^ that the resonance has decayed into 
a stable particle of energy E at time r, 


dP T 

d E 


\{ E \ e - iHT /^ n )\ 2 . 


(2.3) 


As is well known, the survival probability of a resonant (Garnow) state abides by 
the exponential law, 

p s (r) = e -rws . (2.4) 

The initial decay rate p s (0) associated with the survival probability satisfies 


r R = —hp 8 (0) = -h 


dps(r) 


dr 


T —0 


(2.5) 


In analogy to Eq. (12.5(1 . we define the differential decay width of a resonant state 
as follows: 


d T{E) _ dP T=0 _ d /dP T \ 

dE ~ dE dr V dP ) 


( 2 . 6 ) 
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where ^ is the time derivative (rate of change) of The total decay width is 
obtained by integration over the scattering spectrum of the free Hamiltonian, 


r 



dr (E) 
d E 


(2.7) 


This formula is valid when the energy E is the only quantum number needed to describe 
the stable, asymptotic states. If such states are described by additional quantum 
numbers, we just need to sum/integrate over such quantum numbers. 


Three comments are in order here. First, in the standard approach, one uses the 
Golden Rule to calculate the width T^f using Eq. Ill.4)1 . and then one assumes that 
Ti_).f coincides with the width T R that arises from the pole of the 5'-matrix. However, 
as discussed in Appendix 03 the value of the decay width T of Eq. (12.7(1 is in general 
different from the value of T R . Thus, T should be interpreted as another way to 
characterize the strength of the interaction between the resonance and the continuum!]] 


Second, because the probability in Eq. (12.3)1 can be defined also for a square- 
integrable wave function p, we can define the differential and the total decay widths for 
the transition from a state (p into a continuum. We can then expand p in terms of the 
resonant (Garnow) states and a background, | p) = c R |^ R ) + |bg), where |bg) is the back¬ 
ground, and where it has been assumed that the system has only one resonant state 
|z R ). The resulting expression will have a resonant part associated solely with |z R ), 
and other terms that involve |bg). As it usually happens with resonant expansions, the 
term associated exclusively with 1%) carries the contribution from the resonance per 
se. Such resonant contribution is what Eqs. (12.3(1 . (12.6(1 . and (12.7(1 contain. The other 
terms (which, although important, are not the focus of the present paper) contain the 
non-exponential contributions to the decay constant. 

Third, the resonant decay width of Eq. (12.7)1 depends on the normalization of the 
Gamow states. When we use Zeldovich’s normalization, the decay width has dimensions 
of energy, but if we used another normalization, the value and possibly the dimensions 
of T would change. 


Before finishing this section, it may be useful to add a few words about the role 
of resonant expansions in quantum mechanics. In a resonant expansion such as \<p) = 
c R |z R ) + |bg), the resonant state |z R ) is supposed to carry the resonance’s contribution 
to the state p (including the exponential decay), and the background is supposed to 
carry the non-resonant contributions (including deviations from exponential decay). 
When one calculates, for example, the survival probability using a wave function ipi = 
c Rjl |z R ) + |bg x ), one usually obtains an exponential decay for intermediate times and 
deviations from exponential decay at short and long times. However, if we used a wave 
function p 2 = c R , 2 |£ R ) + |bg 2 ) that is only slightly different from p 1: we would obtain a 
slightly different survival probability with slightly different deviations from exponential 


1 Interestingly, in particle physics there are sometimes two and even three definitions of the width 
of unstable particles [1|, and only one of them coincides with the width that arises from the pole of 
the S'-matrix (mg. 
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decay. The question then arises: What is the wave function of the resonance? Is it cpi, 
<P 2 , or something else? Or even worse, do tpi and <p 2 represent the wave functions of 
two different resonances? The Gamow-state description of resonances addresses these 
questions very easily: For each pole of the S-matrix, one can construct a unique Garnow 
state that is identified with the wave function of the resonance. In our example, there 
is only one resonance, and ip\ and (p 2 are two different approximations of one and the 
same resonant state. Thus, by identifying the resonant (Garnow) state with the wave 
function of a resonance, one has a clear way to prescribe what is resonance from what 
is not resonance. 

Identifying the Garnow state with the resonance’s wave function is very similar to 
identifying a monochromatic wave with a plane wave. Both a plane wave and a Garnow 
state are not square-integrable functions, but when you prepare a square integrable 
wave function very close to a plane wave or a Garnow state, you have, for all purposes 
and intends, a monochromatic state or a Garnow state. 


3 Derivation 


Let us assume that our system has been prepared in a resonant state |z R ). Such state 
satisfies the following integral equation mm- 

i 


\Zr) = 


z R — Hn + iO 


V\z R ). 


(3.1) 


By multiplying Eq. (13.ip by e 12RT / fi , and by taking into account that e lz ^ T / h \z n ) = 

~\Ht / 


e ‘- T / h \£ r ), we obtain that 

-i Hr 


e z n )=e 


-i Z R t/H_ 


V\z R ). 


z R — Ho + iO 

The “scalar product” of Eq. (13. 2 p with (E | can be written as 

1 


(E\e- iHT/n \z R ) = e ~ iZRT/h {E\- 


tzV\z r ) . 


When 


i 


z R —Ho+iO 


z R — Ho + iO 

acts to the right, as it does in Eq. (13.31) . we have that 
1 1 


(E |- 


z R - H 0 + iO z R - E 
Substitution of Eq. (13. 4 p into Eq. (j3.3[) yields 


(E |. 


(E\e 


~ iHr / H \ ZR \ — e “ i2 R T / h 


Z R - E 


<£|VK). 


By combining Eqs. (13.5ft and (12.3p . we obtain 


d P T (E) 
d E 


= \(E\e 


-iHT/h\ z \ 12 _ e -r R r/n 


(e - E R y + (r,/2)2 


\(E\V\z R )t 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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The time derivative of this expression reads 


d Pr(E) 
d E 


_ e -r rt/R 


r Jh 


(E - E n y + (T r /2) 2 


\(E\V\zn)f 


(3.7) 


Substitution of Eq. 
state, 


(13.7p into Eq. (12,6[) yields the differential decay width of a resonant 


dr (E) 

d E 


(E - E R y + (r R /2 ) 2 


\(E\V\Zn)\ 2 . 


(3.8) 


Substitution of Eq. (13.8jl into Eq. (12.7j) yields the total decay width of a resonant state, 


r 



{e - e r ) 2 + (r R /2 ) 2 


l<S|r|2a)| 2 . 


(3.9) 


Thus, the decay width of a resonant state depends on the matrix element squared of 
the interaction potential and on the Lorentzian lineshape. The energy domain of the 
Lorentzian, however, does not extend from —oo to +oo, but it is truncated at the lower 
energy threshold E — 0. Equations (13. 8 p and (13. 9 p contain no approximations and, as 
we are going to see below, they bear a strong resemblance with the standard Golden 
Rule of Eqs. (1 1.3 p and (II. 4p in the approximation that the resonance is sharp and far 
away from the threshold. 


In order to relate Eqs. (13. 8 p and (13.9p with the standard Golden Rule, we are going 
to use the well-known result that, when its width tends to zero, the Lorentzian becomes 
the delta function, 


5(x) = lim 


e/7T 


e-»0 X 2 + e 2 


(3.10) 


Let x = (E — E r ) and e = T R /2. Then, when the resonance is sharp (i.e., when 
T r /( 2E r ) is small), and when we can neglect the lower energy threshold (i.e., when the 
resonance is far from the E = 0 threshold), we have that Eq. (I3.10P yields 


S(E - E r ) 


lim -- V ^ 2lT ] , v 

r R ->o (E - E R ) 2 + (T r /2) 2 


By combining Eqs. (13.8ft . (13.91) and (13.111) . we obtain 


d T{E) 
d E 


^\X\V\z R )\ 2 5{E — E r ) , 


(3.11) 


(3.12) 


r«2 I r|(£ ri |l/| ZR )| 2 . 


(3.13) 


Equation (13.12p has the same form as Eq. (1 1.3 p . Equation (I3.13P has the same form as 
Eq. (II.4p . except that in Eq. (13.13ft the “density of states” is 1. 

It is surprising that the decay widths of Eqs. (13.8ft and (13.9ft become very similar to 
the standard Golden Rule in the approximation that the resonance is sharp, since the 
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derivation leading to Eqs. (j3. 12[) and (j3. 12[) is completely different from the derivation 
leading to the standard Golden Rule. Although it is not clear why such is the case, in 
Appendix |B] we explore a possible explanation. 

In the literature, one can find several expressions for the transition probability 
of Eq. (13. 6p and the decay widths of Eqs. (13. 8p and (13. 9p . see Refs. [TTTfT9] . Such 
expressions always involve the matrix element squared of the interaction potential and 
the Lorentzian lineshape. In particular, the decay widths (13.8p and (13.9p coincide with 
those provided in Ref. [T9], although the expressions of Ref. [T9] are derived from a 
time-dependent probability that is different from that in Eq. (I-).01) PI 


4 An illustrative example: The delta potential 


We are going to use Eqs. (13.8p and (13.9p to calculate the decay width of a particle 
trapped by a delta potential [28l 132] located at r = a, 

V ( r ) = gS(r — a ), (4.1) 


where g accounts for the strength of the potential. Since this potential is spherically 
symmetric, we can work in the radial, position representation. For simplicity, we will 
restrict ourselves to s-waves (zero angular momentum), since the higher-order cases are 
treated analogously. The eigenfunctions Xo( r ', E) of the free Hamiltonian Hq satisfy 

-T--T 2 Xo(r-,E) = E X o(r-,E). (4.2) 

2m d r z 

The delta-normalized solution of Eq. (14.21) is well known (see for example Ref. [33]), 


Xo(r]E) 


1 2 m/h 2 


7r 


k 


sin ( kr ) , 


0 < r < oo , 


(4.3) 


where k = 
H satisfy 


subject to 


yjjJrE is the wave number. The eigenfunctions y(r; E) of the Hamiltonian 
h 2 d 2 


2m dr 2 
the boundary conditions 


+ g5(r-a) y(r; E) = Ey(r; E ), 


(4.4) 


x(0; E) = 0, 

(4.5) 

y(a+; E) = y(a—; E ), 

(4.6) 

x'(a+]E) x'(a ; E) — ^_ 2 y(a; E). 

(4.7) 


2 Further variations on Fermi’s Golden Rule can be found in Refs. [20K27] . 
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The solution of Eq. (14.41) subject to Eqs. (14.50 - (14.7|) is given by 


X (r-,E) = A(E) 


sin (At) 0 < r < a , 

Ji(E)e lkr + J 2 (E)e~ lkr a < r < oo , 


(4.8) 


where 


ME) = 


i ka 


Aka 


1 — i , -Jr- ) sin(fca) — icos (ka) 


ME) = — 


1 + i 


h 2 k 

2 mg 
~¥k 


sin (ka) + i cos (ka) 


(4.9) 

(4.10) 


and where A(E) is an overall normalization constant. The 5"-matrix reads 

ME) 


S(E) — 


ME) 


(4.11) 


The resonant states are obtained by imposing the purely outgoing boundary con¬ 
dition on the Schrodinger equation. In our case, this is tantamount to imposing that 

J2(e) = o, 

2777/(7 \ 

' ' sin(fca) + icos (ka) = 0. (4-12) 


1 + i- 


h 2 k 


We will denote the complex solutions of Eq. (14.12)1 by z n , n = 1,2,.... The corre¬ 
sponding resonant states read as follows: 


w(r; Zn) = N n 


sin (k n r) 0 < r < a , 
a < r < oo , 


Jl(Zn) 

k n r 


(4.13) 


where Zeldovich’s normalization constant N n is given by the residue of the S'-matrix 
at the complex resonant wave number k n , N 2 = ires [*S'(q , )] c;= ^• 

The matrix element of the potential is given by 


(E\V\z n ) = / drxo{r]E)gS(r - a)u(r-,z n ) 


= 9\~ 


1 2 m/h 2 


Ti k 


sin (ka) N n e 


i k n a 


(4,14) 


Thus, Eq. (13.8j) yields 

dT n (E) 
d E 


r, 


(E - E n ) 2 + (T n /2) 2 7 rkh 2 


^sm 2 (ka)\N n \ 2 e 2pna 


(4.15) 
























where j3 n is the imaginary part of the complex wave number k n , i.e., k n 
The total decay width is given by 


i fin- 


where 


f _ 2m 9~ \ N \2 2P n a/~i 


° n J 0 dE (E - E n ) 2 + (T n /2) 2 k 


— sin 2 (ka) 


(4.16) 

(4.17) 


Because the delta-shell potential is very simple, one can calculate the decay width 
analytically. However, for more complicated potentials, it may not be possible to 
calculate the resonant states exactly, and one may have to resort to approximations. 


5 Partial decay widths and branching fractions 


When a resonance has different decay modes, one usually defines the branching fractions 
to account for how likely the resonance will decay into each mode. Experimentally the 
branching fractions are measured as the number of particles that decay into a given 
mode divided by the total number of decays. 

In order to obtain the theoretical branching fractions, we need to define the partial 
decay widths. Let us assume for simplicity that the resonance has two decay modes 
and that each decay mode can be described by a single, discrete, quantum number 
j = 1, 2. Similarly to Eq. (j3.6|) . the probability (per unit of energy) that the resonance 
decays at time r into the mode j is 


d P S ,r(E) 

d E 


\{E,j\e~ iHr/h \z R )\ 2 = e- rRT / h 


(E - E n y + (r R /2 ) 2 


I (E,j\V\z n )f 


3 = 1,2. 
(5.1) 


Similarly to Eq. (13.SD , the partial differential decay width of a resonant state is 


dlj(-E) 


r F 


d E (E - £ r ) 2 + (r R / 2) 2 

The partial decay width is 

— f°° T 

Tj = Jo ^ C E - E r ) 2 + (T r / 2) 2 
The resulting branching fractions are 

f, 


\(E : j\V\z r )\ 2 , J = 1,2. 


\(E,j\V\z R )\ 2 , J = 1,2. 


B( R->j) = =-* ir , j = 1,2. 
t 1 + t 2 


(5.2) 

(5.3) 

(5.4) 


Similarly to Eq. (13. 13[) . when the resonance is sharp, the partial decay widths are 
approximated by 

r? « 2 tt|(E r , j\V\z R )\ 2 , J = 1,2. (5.5) 
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j = 1,2, 


(5.6) 


In such an approximation, the branching fractions read 


B (R -> j) « 


1 ( E r , j | H |; r )| 2 

|(£ r , 1 | E | z r )| 2 + |(£ R , 2 | E | z R )| 2 ’ 


which is the same result we would obtain if we used the decay widths of Eq. (11.41) with 
| Ei) = |z R ) and \Ef) = \E R ). Thus, in the approximation that the resonance is sharp, 
the branching fractions obtained from the probability distribution of Eq. (15.1|) are the 
same as those obtained from the Golden Rule. However, when the resonance is not 
sharp, we would obtain different branching fractions. Because the branching fractions 
do not depend on the normalization of the resonant state, Eq. (15.41) will in general yield 
a different result from that of Eq. (15 .6H . no matter what normalization is used for the 
resonant state. This, at least in principle, makes the formalism of the present paper 
falsihable. 


6 The (dimensionless) decay constants 


There are two main kinds of decay experiments. In one kind, one measures the number 
of decay events as a function of time, which corresponds to measuring a time-dependent 
probability distribution such as the survival probability. In another kind, one measures 
the number of decay events as a function of the energy, and one obtains an energy 
distribution of decay events such as the invariant mass distributions of particle physics 
(see for example Refs. [34H42] ). The purpose of this section is to describe the energy 
distributions of decay events in terms of the resonant states. 

Although a scattering experiment is a time-dependent phenomenon, one usually 
describes it using time-independent quantities such as the cross section and the S- 
matrix. Such time-independent quantities can be obtained from solutions of the time- 
independent Schrodinger equation. Similarly, a decay process is a time-dependent 
phenomenon, but it can be described by time-independent quantities such as the energy 
distributions of decay events. The experimental energy distributions of decay events 
are simply the number of decay events per unit of energy as a function of the energy, 


# decay events 
energy 


vs. energy 


( 6 . 1 ) 


When we divide the number of decay events by the total number of events N 0l the 
corresponding experimental distribution corresponds to the theoretical probability dis¬ 
tribution of decay, 


d P(E) / decay probability 

d E \ energy 

where is the (theoretical) probability 

decayed into a stable particle of energy E. 


vs. energy 




1 # decay events 


_ No energy 

per unit of energy that the resonance has 


vs. energy^ , 
( 6 . 2 ) 
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When | E) is delta normalized and |z R ) is normalized according to Zeldovich’s 
normalization, the “scalar product” (E\z R ) has dimensions of (^/energy) 1 . Thus, 
|(-£/|z R )| 2 has dimensions of (energy) -1 , and it is natural to identify it with 


d P(E) 
d E 


\(E\Zn)\ 2 . 


(6.3) 


Obviously, the probability of Eq. (16. 3 p is the time-independent version of the probability 
of Eq. (12.31) . The total probability that the particle has decayed into the continuum is 


P 


r°° d P(E) 

Jo 


d E. 


(6.4) 


Similarly to the definitions of the decay widths dh and T, we can define the differential 
and the total decay constants as follows, 

= (6.5) 

r = P= / 0 ”|(B|z b )| 2 dB. (6.6) 

By using Eq. (13.ip . and by following similar steps to those in Sec. El we can express 
the decay constants of Eqs. (16.51) and (16.61) as 


drg) 

d E 


1 

(E - E R y + (T r /2) 2 


\(E\V\z r )\\ 


r 



i 

(E ~ E R y + (T r /2) 2 


I(^|e|^ r )| 2 . 


(6.7) 

( 6 . 8 ) 


These expressions are very similar to the expressions for the decay widths of Eqs. (13. 8 p 
and (13.9p . In fact, it follows from Eqs. (13.8p . (13.9p . (16.7p and (16. 8 p that 


dr(E) _ l dr(E) 

d E Tr d E 7 


(6.9) 


r - 


r 

r R ’ 


( 6 . 10 ) 


As can be seen from Eq. (16.4p or from Eq. (I6.10p . T is a dimensionless constant that 
can be interpreted as another way to measure the overall strength of the interaction 
between the resonance and the continuum. By contrast, has units of (energy) -1 , 
and hence can be interpreted as the probability density that the resonance decays into 
a stable particle of energy E. 

Similarly to the decay width T, one can define the dimensionless decay constant T 
of a square-integrable wave function ip. When we perform a resonant expansion of <p 
and substitute such expansion in the expression for T, we obtain a contribution that 
is associated exclusively with the resonant state, and other contributions that involve 
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the background. The contribution that is associated exclusively with |z R ), which is 
given by Eqs. (j6.7[) and (16.81) . is the contribution of the resonance per se. It should be 
noted that the (total) decay constant T associated with a normalized wave function p 
is equal to one, as it should be, since the number of events per unit of energy must add 
up to the total number of events. 

Similarly to and Tj, we can also define the partial decay constants as follows: 


dr j(E) 
d E 

r i = 


Pi 


d Pj(E) 
d E 



1 

(E - E K y + (TV2) 2 

1 

(E - E R y + (T r /2)2 


\(E j 3 \V\z r )\ 2 , 

■\(E,j\V\z r )\\ 


3 = 1,2, (6.11) 

i = 1,2, (6.12) 


where Pj is the probability that the resonance decays into mode j. It follows from 
Eqs. (15.2ft . fl5.3p . (16.lip and (16. 12p that 


dVj(E) = l dr j(E) 
d E r R d E ' 



3 = 1 , 2 , 

3 = 1,2. 


(6.13) 

(6.14) 


The main advantage of the partial decay constants Tj over the partial decay widths 
Fj is that the decay constants Tj make the connection between the experimental and 
the theoretical branching fractions more apparent. In order to see this, let Nj be the 
number of decay events into the j mode. Then, the experimental branching fractions 
are defined as 

N- 

s|r ^ jI = m- * = 1 - 2 - (6 ' 15) 

If the experiment is done Nq = N\ + N 2 times (or, equivalently, if there are Nq copies 
of the same resonance), the number of decay events into mode j is equal to the number 
of possible decay events multiplied by the probability that the resonance decays into 
mode j, 

Nj = N 0 Pj = N 0 Tj , J = 1,2, (6.16) 

where we have used Eq. (16.121) in the second step. Substitution of Eq. (16.161) into 
Eq. (16.151) yields 

= 3 = 1,2. (6.17) 

Because of Eq. (16.141) . the theoretical branching fractions of Eq. (16.17P are the same 
as those in Eq. (15.41) . However, in Eq. (15.4p the branching fractions are expressed as 
the quotient of two dimensionful quantities, whereas in Eqs. (16. 15ft and (16.171) they 
are expressed as the quotient of dimensionless quantities. Thus, it seems preferable 
to express the theoretical branching fractions in terms of the (dimensionless) partial 
decay constants Tj rather than in terms of the partial decay widths Fj. 
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7 Conclusions 


The resonant (Gamow) states are the natural wave functions of resonances, and they 
lead to expressions for the differential, partial, and total decay widths of a resonance. 
In this paper, we have presented a simple derivation of the expressions for such decay 
widths and have shown that, in the approximation that the resonance is sharp and far 
from the lower energy threshold, the expressions bear a close resemblance to the stan¬ 
dard Golden Rule. When the potential is simple, one can obtain an exact expression 
for the decay width (as in the example of the delta-shell potential), although for most 
potentials one is bound to use approximations. We have also constructed the branch¬ 
ing fractions for a resonance that has several decay modes. The resulting branching 
fractions coincide with the branching fractions obtained by way of the Golden Rule 
only when the resonance is sharp. Thus, when the resonance is not sharp, the branch¬ 
ing fractions of the present paper can, at least in principle, be distinguished from the 
standard ones based on the Golden Rule. 

We have also introduced the (dimensionless) partial, and total decay constants of 
a resonant state. Using these dimensionless decay constants, we can express both the 
theoretical and experimental branching fractions as the ratio of dimensionless quan¬ 
tities. In addition, the differential decay constant can be identified with the energy 
distributions of decay events, which are a non-relativistic analog of the invariant mass 
distributions of particle physics. 

Although the results of the present paper are restricted to the non-relativistic case, 
it is not hard to imagine a relativistic extension. Essentially, in the usual expression for 
the decay rate [32], instead of a delta function we should have a truncated, relativistic 
Breit-Wigner lineshape. 
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A Relationship between T and T R 


When the survival probability is given by p s (r) = e 
T r are related by 


TrT / h , the initial decay rate and 


Tr = —hp s (0) . 


(A.l) 
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However, if the survival probability was p s (r) = p^e rRT / fi , the initial decay rate and 
r R would be related by 

r R p 0 = -hp s (0). (A.2) 

Thus, only when the probability is given by the exponential function with no prefactor, 
does Eq. (1A.1I) hold. As we are going to see below, because the probability of Eq. (12. 3 I) 
has a prefactor in front of the exponential, the decay width of Eq. (12.71) is different 
from the A-matrix width T R . 

Equation (12.3ft can be written as 


dP, 


d E 


l = e- r ^/«|<£W| 2 . 


Thus, the differential decay width (12.61) can be written as 

d T(E) 


d E 


= r„|{£Wi 2 . 


The total decay width is 


r = r R / dP|(P|^ R )| 2 . 

J 0 

Thus, T would be equal to T R if 

poo 

/ dE \(E\z R )f = 1. 

Jo 

Since 

poo poo 

/ d E \(E\z r )\ 2 = / dr \{r\z, 


2 

R/ I i 


(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 


the decay width T would be equal to the A-matrix width T R if the absolute value 
squared of the Gamow state u{r\ * R ) = (r|z R ) was normalized to one, 


dr |u(r; ^ R )| 2 = 1. 


(A.8) 


However, since Zeldovich’s normalization is (formally) given by 


dr [u{r]z R )} 2 = 


1 , 


(A.9) 


the decay width T is in general different from the A-matrix width T R . 

Finally, we would like to note that, in the standard derivation of the Golden Rule, it 
is assumed that the transition probability per unit of time (decay rate) is independent 
of time. However, the decay rate associated with the the survival probability of a 
resonant state, p s (r ) = — (T R /h)e _rRT//fi , is not independent of time, unless we restrict 
ourselves to a zero-order approximation. The same holds for any other probability 
distribution that follows the exponential decay law. 
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B The analogy with the Golden Rule 


In this Appendix, we would like to point out a possible explanation for why the resonant 
decay widths of Eqs. (13. 8 p and (j3.9[) bear such a close resemblance to the standard 
Golden Rule of Eqs. (11.3p and (jl.4p in the approximation that the resonance is sharp. 

The Lippmann-Schwinger equation for “in” states can be written as 

' £+ > = V+E^kTTe^- ( R1 > 

By writing Go = £ ._^ 0+je , and by successive substitutions, we obtain the Born expan¬ 
sion of the Lippmann-Schwinger equation, 

| E + ) = | E) + G 0 V\E) + (G 0 V) 2 \E) + • • • . (B.2) 

In this expansion, | E) is the term that one would obtain if there was no interaction, 
i.e., if V = 0. Thus, one usually defines the scattered state, \E S ) = \E + ) — | E), as 
the state that contains the actual effect of the interaction. Equation (IB. 2D can then be 
written as 

I E s ) = G 0 V\E) + (G 0 V) 2 \E) + ■■■ . (B.3) 

The first-order approximation to the scattered state is given by, 

\E S ) « G 0 V\E) . (B.4) 

On the other hand, the integral equation satisfied by the resonant states, Eq. (13.ip . 
can be written as 

k> = Go(*r)^I*r> ■ (B.5) 

When the resonance is sharp, one can approximate the resonant state by a bound state 

|£r>, 

\z R ) ^ G 0 (E r )V\E r ) . (B.6) 

The similarity between Eqs. (1B.4D and (1B.6D may be the reason behind the similarity 
between the standard Golden Rule of Eqs. (I1.3p ~ (ll. ID and the approximate decay 
widths of Eqs. (13.12D (13. 13D . 

Although Eqs. (1B.4D and (1B.6D are analogous, there is a major difference between 
them. The approximation in Eq. (1B.4D relies on the assumption that the potential 
V is weak, whereas the approximation in Eq. (1B.6D relies on the assumption that the 
potential V is strong enough that the resonance can be considered nearly a bound 
state. 
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